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ABSTRACT 

The radiation pressure of a radio pulsar does not necessarily disrupt a surrounding disk. The position of the 
inner radius of a thin disk around a neutron star, determined by the balance of stresses, can be estimated by 
comparing the electromagnetic energy density generated by the neutron star as a rotating magnetic dipole in 
vacuum with the kinetic energy density of the disk. Inside the light cylinder, the near zone electromagnetic 
field is essentially the dipole magnetic field, and the inner radius is the conventional Alfven radius. Far outside 
the light cylinder, in the radiation zone, |E | = |B | and the electromagnetic energy density is (S)/c oc 1/r 2 where 
S is the Poynting vector. Shvartsman (1970) argued that a stable equilibrium can not be found in the radiative 
zone because the electromagnetic energy density dominates over the kinetic energy density, with the relative 
strength of the electromagnetic stresses increasing with radius. In order to check whether this is true also 
near the light cylinder, we employ global electromagnetic field solutions for rotating oblique magnetic dipoles 
(Deutsch 1955). Near the light cylinder the electromagnetic energy density increases steeply enough with 
decreasing r to balance the kinetic energy density at a stable equilibrium. The transition from the near zone to 
the radiation zone is broad. The radiation pressure of the pulsar can not disrupt the disk for values of the inner 
radius up to about twice the light cylinder radius if the rotation axis and the magnetic axis are orthogonal. This 
allowed range beyond the light cylinder extends much further for small inclination angles. The mass flow rate 
in quiescent phases of accretion driven millisecond pulsars can occasionally drop to values low enough that the 
inner radius of the disk goes beyond the light cylinder. The possibilities considered here may be relevant for 
the evolution of spun-up X-ray binaries into millisecond pulsars, for some transients, and for the evolution of 
young neutron stars if there is a fallback disk surrounding the neutron star. 

Subject headings: accretion disks — stars: individual SAX J1808. 4-3658, Aquila X-l — stars: neutron — X- 
rays:binaries 



1. INTRODUCTION 

A neutron star can interact with a surrounding disk in a 
variety of modes. These different modes of interaction can 
lead to various astrophysical manifestations and are of key 
importance in classifying neutron stars. The mode of interac- 
tion and hence the evolutionary stage is determined by the 
location of the inner radius of the disk with respect to the 
characteristic radii, the corotation radius R c = (GM/i7 2 ) I//3 
and the light cylinder radius Rl = c/Q, where is the an- 
gular velocity and M is the mass of the neutron star. Start- 
ing with the pioneering work of IShvartsmanl (O]?70a), three 
basic modes of interaction -accretor, propeller and ejector- 
jLipunovll9 92) of a neutron star with a surrounding disk have 
been identified. If the inner radius of the disk is beyond the 
corotation r adius, the system is expected to be in the pro- 
peller stage JShv artsmaiJll970allIl larionov & Sunvaevlll975l 
IDavidson & Ostrikerlll973t lFabiarJll975l) . The eiector (radio 
pulsar) stage is assumed to commence when inner radius is 
also beyond the light cylinder radius. 

The disk should be disrupted in the region where ma gnetic 
and matter stresses are comparable JPring le & Ree ^| | 1972t 
IDavidson & Ostrikerll97"3llLamb. P e thick & Pine sll973l) and 
the inner radius of the disk, R,„, is estimated by balancing the 
kinetic energy density (half the ram pressure) of the disk with 
the electromagnetic energy density (same as electromagnetic 
pressure) of the radiation generated by the neutron star. In or- 

1 Present address: Harvard-Smithsonian Center for Astrophysics, 60 Gar- 
den Street, Cambridge, MA 02138; yeksi@cfa.harvard.edu 
Electronic address: yavuz@sabanciuniv.edu, alpar@sabanciuniv.edu 



der that the equilibrium found in this way is stable, the electro- 
magnetic energy density at r < R,„ should be greater than the 
kinetic energy density, i.e. the electromagnetic energy density 
at Rj„ should increase more stee ply than the ki netic energy 
density in approaching the star (lLipunovl [T992). The usual 
estimate for the kinetic energy density has the dependence 
£k oc r" 5 / 2 (see §3) on the radial distance r from the center 
of the neutron star. In the near zone (R* < r <C Rl) where R„ 
is the radius of the neutron star, electric energy density, £ e , 
is negligible and the kinetic energy density is balanced essen- 
tially by the magnetic energy density, £ m , of the dipole field 
of the neutron star at the Alfven radius. This is a stable equi- 
librium point because £ em oc r" 6 is steep enough to balance 
Ek oc r" 5 / 2 i.e. £ em > £ K for r < /?,„ and £ K > £ em for r > Rj„. 
In the radiation zone (r ^> Rl) where electromagnetic energy 
density £ em oc r" 2 , a stable inner boundary can not be found 
because then £k would be greater than £ em for r < and 
£em > £k for r > Rj„, and the d isk is disrupted. This was first 
noticed by Shvartsman ( 1970b) who with the above argument 
concluded that the stable equilibrium outside the light cylin- 
der wo uld only be po ssible beyond the gravitational capture 
radius (Lipuno vl 19921) . 

Usual estimates employ a rotating magnetic dipole in vac- 
uum to generate the electromagnetic fields of the neutron star. 
In the near zone (/?* < r <C Rl) the magnetostatic energy den- 
sity B 2 /8ir is employed, while in the radiation zone (r ^> Rl), 
at distances far outside the light cylinder, the radiation pres- 
sure (S) /c, where S is the Poynting vector, is used. For sim- 
plicity one usually employs a piecewise defined electromag- 
netic energy density which scales as r" 6 for r < Rl and as r~ 2 
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FIG. 1 . — Two possible equilibrium configurations around a rotating dipole. 
Equilibrium at the near zone (point A = Alfven radius) is stable because mag- 
netic energy density drops more rapidly than the kinetic energy density going 
beyond this point. Equilibrium at the radiation zone (point B) is not stable 
because the magnetic pressure drops less rapidly than kinetic energy density 
and disrupts the disk beyond this point. With a piecewise defined field con- 
figuration as above, one has to conclude that the disk has to be disrupted once 
the inner radius goes beyond R^. (This figure is not included in the accepted 
paper. It is added for explanatory purposes. ) 
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FIG. 2. — Radial structure of the electromagnetic energy density around 
the pulsar. Solid (colored red in electronic edition) line is the electromag- 
netic energy density, long dashed (colored green in electronic edition) line is 
the electrical energy density E 2 /8n, short dashed (colored dark blue in elec- 
tronic edition) line is the magnetic energy density B 2 /87r, and dashed and 
dotted (colored light blue in electronic edition) line is (S r )/c. Magnetic en- 
ergy density dominates inside the light cylinder radius. Far outside the light 
cylinder electric and magnetic energy densities become equal and {S r )/c is 
equal to their sum. See the electronic edition of the Journal for a color version 
of this figure. 



for r > Rl- With such a model one necessarily concludes that 
the disk will be disrupted if the inner radius of the disk goes 
beyond the light cylinder radius (see Figure^. 

In this work we present the full expressions for electromag- 
netic pressure around a rotating dipole in vacuum, reducing 
to the conventional expressions in the near zone and the ra- 
diation zone, and describing the transition across the light 
cylinder. We use these results to investigate the location and 
stability of the Alfven radius and inner radius of thin Keple- 
rian disks. This information is relevant if the mass flow rate 
through the disk declines or the neutron star rotation period 
evolves such that the inner radius is beyond the light cylinder. 
According to the convent ional pulsar magnetophere models 
JGoldreich & JuliarT l969). the pulsar activity can commence 
if the inner radius is outside the light cylinder. According to 
the simple piecewise model described above, the disk would 
be swept away. In the present model with a rather broad tran- 
sition from the near zone to the radiation zone, there will be a 
wide range of allowed values for the mass flow rate for which 
the inner radius of the disk is outside the light cylinder while 
the disk can survive the radiation pressure. In this case, the 
pulsar activity will commence but depending on the inclina- 
tion angle, the presence of the disk may destroy the coherent 
pulsed emission in the radio band. 

In its quiescent stage, the luminosity of the accre- 
tion driven millisecon d pulsar SAX J1808.4-3658 

( Wiina nds & van der Klisl 119981) drops below L x ~ 10 32 
erg s" 1 dCampana et al J 12002). Assuming the system passes 
through a propeller stage the mass inflow rate, which deter- 
mines the inner radius, can be much greater than the mass 
accretion rate on to the star determining the luminosity. As 
we do not know what fraction of the inflowing mass can 
accrete in the propeller stage, we can not estimate the inner 
radius from the observed luminosity. Assuming 1 per cent of 
the inflowing mass accretes on to the star, we would estimate 
an inner radius for the disk greater than the light cylinder. 



iBurderi et al.l J2003) suggested that the optical properties of 
this sour ce in the quiescent s tage indicate an active pulsar 
(see also Campanaetal. (2004)). We find that such a disk 
can survive the radiation pressure of a turned on rotation- 
powered pulsar . This may be the case in some stag es of 
transients (e.g. Camp ana et alJ IT998): Zam anovl ( 1995)) and 
of the evolution of young neut r on stars if they have fa llback 
disks JMichel & Desslerlll98U iMarsden et alJl20Qirl 120021: 
iMenou et alJ200UlAlpar et alJ200UlEksi & AlpaJ2003l) . 

In §2 we derive the electromagnetic ener gy density from 
the global electromagnetic field solution of iDeutschl ( 119551) 
for obliquely rotating magnetic dipoles. In §3 we derive the 
inner radius of the disk. In §4 we discuss implications for the 
accretion driven millisecond pulsar SAX J 1808 .4-3658 which 
is probably representative of the millisecond X-ray pulsars, 
for certain other transients and for fallback disks around radio 
pulsars. In §5 we discuss our results. 

2. ELECTROMAGNETIC ENERGY DENSITY 

A rotating dipole can radiate electromagnetic energy if the 
spin and magnetic axes are not aligned. A global solution for a 
perfectly conducti ng, rig idly rotating star in vacuum was first 
given by Deutsch ( 1955). The relevanc e of this solution fo r 
radio pulsars is thoroughly dis cussed by [ Michel & Li ( 1999). 

The early papers on pulsars JOstriker & Gunnll9 68: Pacini 
[T967| U%8F iMichel & Goldwird Il970t iDavis & Go ldstein 
1970) employed the vacuum solutions of Deutsch (1955). 
The vacuum assumption can not actually be realized for neu- 
tron stars because the magnetospheric field itself is strong 
enough to tear off electrons fro m the surfac e of the neu- 
tron star. Indeed, iGoldreich & Juhanl (f!%9T) showed that 
an axisymmetric pulsar magnetosphere should have a coro- 
tating plasma described by the hydromagnetic approxima- 
tion E + v x B = with a corresponding charge density V • E. 
Whether the pulsar electrodynamics is dominated by the 
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iGoldreich & Juliar] i ll 9691) currents, or by the vacuu m fields 
modified by rather smaller currents is still of debate ( Michel 
119911) and a global s elf-consistent solu tion of the problem is 
yet to be presented. Kaburaki (1981) argued that the exis- 
tence of a corotating plasma does not alter the field struc- 
ture drastically because it simply means an effecti ve incr ease 
in the radius of the conducting star (see Melatos (119971) for 
an application of this idea) using the Deutsch ( 1955) solu- 
tion. The presence of the di sk might alter the field structure 
(e.g. Eisner & Lamb ( 1977)) even more than the corotating 
plasma and further complicat e the magnetospheric currents 
( Treves & Szuszkiewics 1995). In the absence of a consis- 
tent solution of the problem of the pulsar magnetosphere with 
corotating plasma w i th or without a disk, we, for simplicity, 
employ the lDeutschl ( 119551) solutions for a rotating dipole in 
vacuum. We quote the fie ld solutions of Deutsch] ( 119551) . as 
corrected by Michel & Li] ( 119991) . in the Appendix, and use 
them in our estimation of the electromagnetic energy den- 
sity at the disk plane. As Maxwell's equations do not pre- 
serve their form in rotating coordinates the problem of rotat- 
ing dipoles is intrinsically relativistic. These solutions can 
be shown to be the limiting case for slow rotation (see e.g. 
Beli nskietalJlfT99l l. 

As obtained in the Appendix, the radial component of the 
Poynting vector S = {c/4tt)E x B at the disk plane (zenith an- 
gle 8 = 7r/2) averaged over a stellar period is 



c 



sin 2 £ 



8nRt 



(D 



Here £ is the angle between the magnetic moment fi and the 
rotation rate ft of the neutron star, and x = r/Ri- Note that 
Eq. ([0 is accurate only for a = R*/Rl <C 1. For the fastest 
millisecond pulsar with P = 1.5 ms, a = 0.14, thus a« 1 
holds for pulsars; in this case (S r ) jc oc r~ 2 for r greater than a 
few R*, (see Eqs. jAl3l -( lAT3t ). 

The dynamical disk time-scale is much longer than the pe- 
riod of the neutron star. The disk would only "see" the average 
field of many stellar rotations. We average the squared fields 
over one stellar period P as (B 2 ) = (1 /P) J* Bjdt where 5; is 
any component of the magnetic field, and obtain the magnetic 
energy density £ m = (fi 2 ) /8 7r as (see Appendix) 



f 



87ri$ 



cos 2 £+-(x 4 + 3x 2 + 5) sin 2 £ 



(2) 



for a -c 1 ■ 

Similarly, in the a <C 1 approximation, the electric energy 
density £ e = (£ 2 )/8 7r is (see Appendix) 



Se- 



ll 



- cos 2 £+i (x 2 + l) sin 2 £ 



(3) 



The total electromagnetic energy density can be found by 
summing up the magnetic and electric energy densities. From 
equations and 0, one finds 



SnRl 



l + -x cos £ + 



x* + 2x z + - ) sin^ 



(4) 

Note that this reduces to £ em oc r 6 for r <C Rl and £ em oc r 2 
for r ^> R L , as expected. We show £ em together with £ m , £ e , 
and (S r /c), all scaled with f^/SirR^, in Figure (|2}- If the sole 
effect of the corotating plasma is to increase the radius of the 
region beyond which vacuum starts from R* to a larger value 




FIG. 3. — Power-law index of the electromagnetic energy density for a 
variety of inclination angles. The power-law index changes between -6 and 
-2. The transition becomes broader for small inclination angles. See the 
electronic edition of the Journal for a color version of this figure. 



R v (lMelatos|[T 997). then a = R v /Ri can have values close to 
unity. In this case the expression for the total energy density is 
more complex than equation @, the coefficients having terms 
with powers of a. The behavior of £ em and its constituents do 
not change qualitatively for a < 1 if we retain such terms. 

The local power-law index of electromagnetic energy den- 
sity, 7 = dln£ em /dlnx, calculated from equation (|4j is 



7 = -6 + 



(8 /9)x 2 cos 2 £ + 4x 2 (x 2 + 1 ) sin 2 £ 
[ 1 + (4/9)x 2 ] cos 2 £ + (x 4 + 2x 2 + 5/2) sin 2 £ 



(5) 



We show 7 for a variety of £ in Figure (|3j- 
The value of 7 at the light cylinder, 7£ 
equation (jSJi, is 



7(x =1), from 



7l =-6 + - 



(8/9)cos 2 £ + 8sin 2 £ 



(6) 



(13/9)cos 2 £+(ll/2)sin 2 £ 

and is plotted in Figure (0J. The kinetic energy density scales 
withx" 5 / 2 (see the following section). Note that for all inclina- 
tion angles jl < -5/2 i.e. the electromagnetic energy density 
is steep enough to balance the kinetic energy density of the 
disk at the light cylinder. 

The radius beyond which 7 > -5/2 can be found by solving 
equation (0 for 7 = -5 /2. We find this to be 



V6 
6tan£ 



4+18tan 2 £+ ^16 + 396 tan 2 £+954 tan 4 £. 



(7) 

Figure[5]shows that the disk inner radius will be stable beyond 
the light cylinder, for distances up to many Rl, practically for 
all radii for aligned (£ = 0°) rotators and up to r„„ = 3Rl for 
orthogonal (£ = 90°) rotators. 

3. THE INNER RADIUS NEAR THE LIGHT CYLINDER 

We estimate the inner radius as the electromagnetic radius 
R em , determined by £ em = £k taking the conventional estimate 
£k = jP\ v \ v ff where p\v\ =M /Anr 2 is from the conservation 
of mass in spherical coordinates, M being the mass flow rate, 
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FIG. 4. — Power-law index of the electromagnetic energy density at the 
light cylinder (x =1), changing with the inclination angle. 
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FIG. 5. — The critical radius, in terms of light cylinder radius at which 
7 = —5/2, as given in equation 



and Vff = y/2GM/r is the free-fall velocity. This gives 



e MV2GM _ 5/2 

OK = TTT- X 



8ttR 



5/2 



(8) 



The electromagnetic radius can be found using equations 10} 
and©: 



X ? / 2 =r 7/2 



l + -x 2 

1 T g A em 



cos 2 ^+[x 4 em + 2x 2 em +- )sin 2 e 



Here xa = Ra/Rl and 

R A : 



(9) 



\mV2GM 



2/7 



(10) 



is the Alfven radius. 

The solution of equation (|9j for x em (xA) is given in Fig- 
ure (|6) for a variety of £. It is not possible to find a solution 
for every xa, corresponding to the fact that at low mass trans- 
fer rates, the kinetic energy density of the disk can not match 



FIG. 6. — The variation of x em with Xa as determined from Eq. |9j, for a 
variety of inclination angles. See the electronic edition of the Journal for a 
color version of this figure. 



the electromagnetic energy density at any radius, and the disk 
will be ejected. The lowest mass flow rate for which the disk 
will not be ejected corresponds to the maximum x A {x em ) as 
seen in Figure 0. From dx A ldx em = follows the solution 
as expected. Now using this back in equation l|9} one finds 
the critical Alfven radius 



X A .crit — Xc 



( 1 + l x ii<) c °s 2 £ + (4rit + 2x ii, + |) sin2 6 



-2/7 



corresponding to the single solution x e 
ical mass flow rate 



(11) 

= x cr i t and to the crit- 



-7/2 



V2GMRl /2 



(12) 



For higher mass flow rates there will be two solutions for x em 
the smaller of which is the stable one. This stable solution 
is beyond the light cylinder for a range of mass flow rates 
M cr u <M < Ml where Ml is the flow rate for which the inner 
radius will be at the light cylinder, x em = 1 . This will happen 
when Xa attains the value 

/13 11 \ -2 / 7 

x a ,l= f yCos 2 ^+ysin 2 n (13) 

as can be obtained from equation (|9jl- Hence 

,2 



M L 



-1/2 



V2GMR 



7/2- 



(14) 



is obtained. Note that both M„„ and Ml are strongly depen- 
dent on the period and the magnetic moment of the neutron 
star (Ml oc M cn - ( oc P" 7 / 2 // 2 ) and their ratio is only dependent 
on the inclination angle £. The range of mass flow rates and 
inner disk radii for which the disk is beyond the light cylinder, 
and stable, is smallest for £ = 90°; for which R cr u = 3Rl and 
M crit = 0AM L . 

Within the model the disk will survive beyond the light 
cylinder, even if the radio pulsar activity turns on. As long 
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as M remains larger than M C nt (i-e. Ra < RA,crit), the disk will 
not be ejected, and the radio pulsar may turn off again when M 
increases to M L , switching back to the propeller phase. The 
threshold for quenching the radio pulsar will reflect loading 
of all closed field lines that lead to the inner gap with plasma 
from the disk. This does not necessarily happen at R em = Rl- 
If the transition to radio pulsar activity happens at some radius 
R t (0 < Rl, then the transient radio pulsar phase with the disk 
not ejected prevails for M < M,, corresponding to R em > R t . 

The scaling we have used which reduces to the Alfven ra- 
dius for r <C Rr presumes a spherical accretion flow. The 
IShakura & Sunvaevl (119731) solutions for thin Keplerian disks 
asymptotically imply that the density of the gas in the disk 
scales as p oc r~ 15//8 . Hence Ek oc pv\ oc pr~ x oc r~ 23 / 8 . This 
scaling (7 = -23/8 = -2.875) is slightly more steep than what 
we have used 7 = -5/2. For a Shakura-Sunyaev thin disk, our 
estimate for the critical inner radius will be somewhat smaller 
than the estimate given above. From Eq. |7} the minimum 
critical radius (£ = 90°) is x c „>, m ,„ = 2.85 while for 7 = -2.875 
we would obtain x cr u m i n = 2.13. 

4. APPLICATIONS 

4. 1 . SAX Jl 808.4-3658, Other Millisecond X-ray Pulsars 
and Transients 

The rotation pe riods of accretion driven millisecond pulsars 
(Wiinands 2004) ra nge between 2.3 ms to 5.4 ms . We take 
the first discovered l lWiinands & van der Klisll998l) accretion 
driven millisecond pulsar SAX J1808.4-3658 with a rotation 
period P = 2.5 ms as an example of this class. The upper 
limit on the magnetic field of the neutron star is 5 x 10 8 Gauss 
JPi Salvo & Burderil2003l) . We assume the magnetic moment 
to be a = 3 x 10 26 Gauss cm 3 . Using these parameters in equa- 
tions O and <E4} we obtain the results shown in Figure0for 
SAX J1808.4-3658. For mass flow rates greater than 10 13 g 
s" 1 for an inclination angle £ = 10° and greater than 2 x 10 14 
g s" 1 for an inclination angle £ = 90°, the source can sustain 
a disk beyond the light cylinder and thus can be in a rotation 
powered pulsar phase. The region between the two curves 
represents the mass inflow rate for which the electromagnetic 
radius is larger than the light cylinder while the disk can sur- 
vive the pulsar activity. At mass inflow rates Ml greater than 
a few 10 14 g s" 1 for all inclination angles, the disk inner radius 
will protrude into the light cylinder, and magnetosperic pulsar 
activity will be quenched by the presence of the disk to the ex- 
tent that all field lines are loaded with plasma from the disk. 
The source will probably be in a propeller phase with only a 
fraction of the mass inflow possibly being accreted while the 
rest is diverted into an outflow. When the inner radius of the 
disk is of the order of the corotation radius, accretion of the 
full mass flow rate is expected. We estimate mass inflow rates 
of the order of 10 16 g s" 1 for this transition for all inclination 
angles. The maximum luminosity of SAX J1808.4-3658 in 
outburst is ~ 2 x 10 35 erg s , which corresponds to a mass ac- 
cretion rate of 10 14 g s" 1 . The discrepancy is probably related 
to our scaling of the disk inner radius with the Alfven radius, 
and the uncertainties in the determination of the Alfven radius 
as well as to the simplifications of the present model in terms 
of the vacuum dipole solutions. In quiescence the luminosity 
drops to ~ 10 32 erg s" 1 (Cam pana et al.ll2002l) . It is possible 
that the mass inflow rate in the disk is large enough that the 
inner radius of the disk is inside the light cylinder but only 
a very small fraction (< 0.1 percent) of this inflowing mass 
accretes onto the neutron star, producing the observed X-ray 



luminosity. iBurderi et alJ (120031) argued that the irradiation of 
the companion by the switched on magneto-dipole rotator, in 
the quiescent stage, can explain the modula tion of the flux in 
the optical (see also Campana et al. (2004)). If the source is 
in the propeller stage accretion of only a small fraction of the 
disk inflow through a limited bunch of field lines will proba- 
bly allow the magnetospheric gaps and pulsar radiation to sur- 
vive. This would be easier for the outer gaps. Magnetospheric 
voltages in millisecond pulsars are of the order of those in the 
Vela pulsar, so optical, X and gamma ray pulsar activity may 
be possible if the outer gap can survive. For small inclination 
angles accretion near the magnetic polar caps would also be 
near the rotational pole and therefore avoiding the centrifugal 
barrier. For paths avoiding the inner gap even radio pulsar ac- 
tivity might survive in the propeller phase, with the disk pro- 
truding inside the light cylinder. Depending on the beaming 
geometry searches for radio and high energy pulsar activity 
might yield very interesting results. SAX J1808.4-3658 was 
selected for this discussion as it is the millisecond X-ray pul- 
sar with the most detailed information (Wiinands 2004). This 
discussion applies to the other millisecond X-ray pulsars also. 
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FIG. 7.— The mass flow rate range for SAX J1808.4-3658 for which the in- 
ner radius is outside the light cylinder while the disk is not swept away by the 
radiation pressure. The solid line corresponds to Ml (eqn. 1141 and the dashed 
line corresponds to M cr it (ecm. 1121 . The domain between the curves corre- 
sponds to the allowed mass flow range for which the inner radius is outside 
the light cylinder while the disk is not disrupted. The observed luminosity 
may be much lower than that would correspond to these mass flow rates be- 
cause the centrifugal barrier will not allow all inflowing mass to accrete onto 
the neutron star. For small inclination angles the allowed mass flow rate range 
increases. 

Other transient low mass X-ray binaries (LMXBs) may also 
be in phases that include excursions of the inner disk towards 
and beyond the light cylinder and back. As an example, the 
luminosity of Soft X- ray Transient Aquila X-l was observed 
Campana et al. ( 1998) to decay from ~ 10 36 erg s" 1 to ~ 10 33 
erg s , from outburst to the quiescent stage. As the system in 
the quiescent stage is likely to be in the propeller mode, the 
observed luminosity of ~ 10 33 erg s" 1 reflects a small fraction 
of the inflowing mass that accretes. The present discussion 
adds to the possibilities that of occasional rotation powered 
pulsar activity. The uncertainty in rotation period and mag- 
netic field does not allow a more detailed discussion in terms 
of our model. Further, the majority of LMXBs do not exhibit 
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rotation periods in their accretion and propeller phases for rea- 
sons of selection in the source, and this is probably also the 
case for most transients. 

4.2. Implications for Fallback Disks 

Neutron stars are born in supernova explosions. It is a 
possibility that some of the eje cted matter fails to a chieve 
the es cape velocity and falls back ( Wooslev, Hege r & Weaver I 
120021) . As the progenitor is rotating, the fallback matter 
might have enough a ngular momentum to settle into a disk 
dMineshiee et alJU997|). The idea of f allback disks around 
pulsars dat es back tolMich el & Desslerl (11981 1). The fallback 
disk m o del JChatteriee et alJ200(]HAlparl200UlMarsden et alJ 
2001a: | E ksi <fe Alpar| |2QQ ^|) of a nomalous X-ray pulsars 
( Mereghet ti et aljr2002t lKaspi 2004) re newed interest in pos- 
sible implicati ons for radio pulsars llMichel & Dessle nll981t 
Marsd en et ail 12001b! 120021: iMenou et al] 120011 Ulnar et alJ 
1200 11) . More recently Black man & Pernal <2004l) suggested 
that the observed jets of Crab and Vela may be collimated by 
fallback disks. 

Young neutron stars might follow a variety of evolution- 
ary paths depending on the initial mass of their fallback disks 
( Alpar 120011) . The evolutionary stages are determined by the 
mode of interaction discussed in § 1 . As the mass of a fallback 
disk is not replenished, mass flow rate in the disk declines 
and the inner radius of the disk, as a consequence, moves 
out. As the neutron star is initially spinning down with pro- 
peller torques, the light cylinder is also moving out, though 
less rapidly than the inner radius. At some stage the inner 
r adius will catch up w ith th e light cylind er. 

IMenou et al. (200l]) and lAlpar et alJ d2001l) assumed that 
the inner radius of fallback disks tracks the light cylinder ra- 
dius. In such a model the disk can assist the magnetic dipole 
radiation torque without quenching the radio pulsar mecha- 
nism. These models would n ot wo rk (Li 2002) if the disrup- 
tion of the disk takes place (Shvartsman 1970b) as soon as 
its inner edge reaches the light cylinder from within the near 
zone. The present work shows that such a disk need not be 
necessarily disrupted till its inner radius moves to a few times 
the light cylinder radius. Whether the inner radius can ac- 
tually track the light cylinder for a long epoch, rather than 
moving beyond the light cylinder into the radiation zone, will 
be studied in a subsequent paper. 

5. CONCLUSION 

We have shown that the electromagnetic energy density of 
a rotating dipole makes a rather broad transition across the 
light cylinder from the near zone dipole magnetic field to the 



radiation zone, the transition being broader for small inclina- 
tion angles, £. For power-law indices 7 for the r dependence 
of the electromagnetic energy density, 7 > -5/2, the regime 
in which the radiation pressure would eject the disk, is not 
attained until r exceeds a few Ri depending on the inclina- 
tion angle between the spin and magnetic axes (see eqn. Q 
and Figure [SJi. Disks that start their evolution with the inner 
radius far inside the light cylinder may evolve to the situa- 
tion where the inner disk radius reaches the light cylinder. In 
this case and even further the electromagnetic stresses of the 
fields generated by the neutron star can balance the material 
stresses from the disk. Near the light cylinder the power-law 
index of the electromagnetic energy density ji varies between 
-5.4 (for £ = 0°) and -4.5 (for £ = 90°) (see Eq. © and Fig©. 
Assuming a thin Shakura-Sunyaev disk with Ek oc r" 23 / 8 does 
not change these conclusions qualitatively and a stable inner 
radius can be found at least up to = 2 times the light cylinder 
radius for orthogonal rotators while this critical radius inside 
which stable equilibrium can be found may extend up to a few 
ten times the light cylinder for small inclination angles. 

When the inner radius of the disk is outside the light cylin- 
der, the pulsar activity is likely to turn on. The presence of 
the disk may effect the coherent radio emission. The mag- 
netic dipole radiation torque will act on the neutron star 
even if a radio pulsar is not observed. Sources on their 
evolutionary path to higher or lower mass inflow rates will 
make a transition from or to a stage with a disk inner ra- 
dius stablely placed in the radiation zone, and with possi- 
ble rotation powered activity. Millisecond X-ray pulsars are 
likely examples of late stages in the evolution of LMXBs 
into millisecond radio pulsars through spin-up by ac cretion 
( Alpar et al1ll982t lRadhakrishnan & Srinivasan 1982). Such 
sources might hover around the transition, exhibiting transient 
behaviour. The stable presence of the disk outside the light 
cylinder for a wide range of mass inflow rates makes repeated 
transitions and sustained transient behaviour possible not only 
around the transition between propeller and accretion phases, 
but also around the transition of the inner radius of the disk 
across the light cylinder. 
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APPENDIX 

THE DEUTSCH (1955) SOLUTIONS 

For an inclination angle £, equatoria l magn etic field strength Bq, and stellar radius R*, the solutions in spherical coordinates 
(r, 9, 4>) are the following ( Michel & Li ( 1999]), Eqns. 89 & 90): 

^3 

B r = 2Z?o — f {cos £ cos 9 + sin £ sin 9 [d\ cos ijj + di sin ip] } , 
Be = Bo—!f {cos£sin#- sin £ cos 6* [(q\ +c/3)cosf/; + (q , 2 + ^4)sin'0] } . 
B^ = Bq—j- sin£ {- \qz cos 29 + d^\ costjj+ [qi cos2(9 + fl?3] sin-0} 



(Al) 



and 



c x L 



Ef> = 



Eq a 2 



-cos£ + — cos £(1 — 3 cos - 9) — j sin£ sm29[q\ cosip + q2 sin?/;] 

2 



cos £ sin 29 + sin £ [(^3 cos 29 - d\ ) cos tp + (qn cos 29 - di) sin tp] 



■ j sm £ cos ^ [(94 - di) cos %j) - (q3 - d\ ) sin tp\ 



Here 



ip = cj)—£lt+x — a 



where 51 is the angular velocity of the star, 
and 

c 

For usual rotation rates a< 1. Moreover, 
and the non-constant coefficients are 



Rl 



R* 
R~l 



R*n 



: 6.28 x 10" 3 



100ms 



Eo = £lR*Bo = caBo 
ax+l , x — a 



d,= 



a 2 +l' 

1 +OX — X 2 



a 2 +l' 

(x 2 — l)a+x 

c/4 = 



and 



01 = - 

02 = 
03 

04= " 



a 2 +l ' " a 2 +l 

3x(6a 3 - a 5 ) + (3 -x 2 )(6a 2 - 3a 4 ) 
a 6 -3a 4 + 36 ' 
(3 - x 2 )(a 5 - 6a 3 ) + 3x(6a 2 - 3a 4 ) 
a 5 -3a 4 + 36 ' 
(x 3 - 6x)(a 5 - 6a 3 ) + (6 - 3x 2 )(6a 2 - 3a 4 ) 
x 2 (a 6 -3a 4 + 36) ' 
(6 - 3x 2 )(a 5 - 6a 3 ) + (6x -x 3 )(6a 2 - 3a 4 ) 



(A2) 

(A3) 
(A4) 

(A5) 
(A6) 

(A7) 



x 2 (a 6 -3a 4 + 36) 



(A8) 
(A9) 
(A10) 



8 EK§1 & ALPAR 

We assume the disk to be at the 8 = ir/2 plane where the field components given in JA1I and iA2\ simplify to 

B r = - sin £ (di cos ip+d 2 sin ip) , 
R l x 



(All) 



B 4> = T^~T sin £ [(?2 - ^4) cos V> + (f^3 - <7 1 ) sin ^1 
,k?x j 



where fx = Z?o/?J and using ( IA6> in (IA2> 



. 2 a 2 



R\x 2 



R\x 2 



sin £ [(^3 + d\ ) cos ^ + (^4 + c/2) sm V'] 1 



(A12) 



£<a = 0. 



We ob tain the radial component of the Poynting vector S = (c/47r)E x B at the disk plane (9 = tt/2) from equations JA1 1> and 
( IA12I . Averaging S r = S • f over a period, we obtain 



S r (i 2 sin 2 £ 



s(x) 



where 

This can be approximated as 



v c ' 8ttR 6 l x 5 
sO) = -[(<73 + di)(q 2 - d 4 ) + (q4 + d 2 ) (d 3 -qi)] 



s(x) : 



"3?" 



forx ~ a; 

for x > a few a. 



(A13) 
(A14) 

(A15) 



We average the square of the fields over one stellar period P as (G 2 ) = (1 /P) J Q Gjdt where G, is any component of the electric 
or magnetic field. 



f - . 



(B 2 



where 



^(cos 2 e + I/ lW sin 2 £ 



fi (x) = (q 2 -d 4 ) 2 + (d i -q\) 2 + Ad 2 l + 4dj 



In the limit a <C 1, this simplifies as /i(x) = x 4 + 3x 2 + 5 and we obtain 



m ~ 87ri$ 



cos 2 £+-(x 4 + 3x 2 + 5)sin 2 £ 



Similarly the electric energy density is 



8- 



8tt/?£ 



2 a 2 Y 1 



^ ( 3 + ^2" ) +^/2Wsin z £ 



where 

and in the limit a<l this simplifies as /2(x) = x 2 + 1 and we obtain 



/2OO = (q3+di) 2 + (q 4 +d 2 ) 2 



-cos 2 £+- (x 2 +l)sin 2 £ 



(A16) 
(A17) 

(A18) 
(A19) 

(A20) 
(A21) 

(A22) 
(A23) 



